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Abstract. Our purpose in this paper is to analyze the Pais-Uhlenbeck (PU) oscillator using 
complex canonical transformations. We show that starting from a Lagrangian approach we 
obtain a transformation that makes the extended PU oscillator, with unequal frequencies, 
to be equivalent to two standard second order oscillators which have the original number 
of degrees of freedom. Such extension is provided by adding a total time derivative to the 
PU Lagrangian together with a complexification of the original variables further subjected 
to reality conditions in order to maintain the required number of degrees of freedom. The 
analysis is accomplished at both the classical and quantum levels. Remarkably, at the 
quantum level the negative norm states are eliminated, as well as the problems of unbounded 
below energy and non-unitary time evolution. We illustrate the idea of our approach by 
eliminating the negative norm states in a complex oscillator. Next, we extend the procedure 
to the Pais-Uhlenbeck oscillator. The corresponding quantum propagators are calculated 
using Schwinger's quantum action principle. We also discuss the equal frequency case at the 
classical level. 
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1 Introduction 



Systems with higher order time derivatives (HOTD) have been studied with increasing interest 
because they appear in many important physical problems. In first place these systems were 
considered to improve the divergent ultraviolet behavior of some quantum field theories [T]. 
However, their energy turned out not bounded from below [2], then involving ghosts and making 
the theory non-unitary [3]. In spite of such drawbacks, higher order derivative theories were 
studied to learn on their improved renormalization properties. Even a renormalizable higher 
order quantum gravity theory was advanced along these lines [4] and the unitarity of a lattice 
form was studied in [5]. See also O El El EJ [lOl E] and references therein for other examples. 
As for dealing with ghosts, attempts can be classified according to whether the approach is 
perturbative [l2l [131 E! or not [El [H [8], but a definite answer is yet to be found. 

*This paper is a contribution to the Proceedings of the Vllth Workshop "Quantum Physics with Non- 
Hermitian Operators" (June 29 - July 11, 2008, Benasque, Spain). The full collection is available at 
|http: //www.emis.de /journals/SIGMA/PHHQP2008.html , 
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Since single particle quantum mechanics can be seen, in the free field or weak coupling 
limit, as a mini-superspace sector of quantum field theory where the spatial degrees of freedom 
have been frozen, it is suggestive to test new ideas on HOTD theories, by using quantum 
mechanics as a laboratory. The model example to do so is the Pais-Uhlenbeck (PU) oscillator |2|, 
which consists of a one dimensional harmonic oscillator Lagrangian plus a term quadratic in 
acceleration. A standard quantum treatment of its degrees of freedom gives rise to a spectrum 
not bounded from below since it consists of the difference of two quantum harmonic oscillators 
spectra [2]. An interesting solution to such difficulty has been recently proposed in [TBI [TTl [TS] 
(see also \19\ [20l \2T\ [22] for a different approach) by finding a quantum transformation which 
gives its Hamiltonian a non-Hermitian, VT symmetric form. In this approach a technique has 
been developed to obtain a physical inner product and thus a Hilbert space [23j . 

The aim of this paper is to show that it is possible to tackle the problem of quantizing an 
extension of the PU oscillator within a Lagrangian and a canonical ormulation, using complex 
canonical transformations and non-Hermitian variables together with reality conditions. This 
approach is motivated by a previous proposal to build complex canonical variables for non- 
perturbative quantum canonical general relativity [241 125] . which achieved the nontrivial task of 
making polynomial the constraints of general relativity, whereas the so called reality conditions 
to be fulfilled by the non-Hermitian variables were proposed to determine the inner product of 
Hilbert space. For canonical general relativity such reality conditions imply the metric of space 
and its rate of change are real quantities. Some effort was devoted to exploit these complex 
variables (see for instance [261 [27] 128] ) however the use of real variables allowed finally for 
important progress [29] [30]. 

The paper is organized as follows. In Section [2] we review the essential aspects of the general 
construction with special emphasis on the boundary conditions. In Section [3] we use as a toy 
model a complexified harmonic oscillator [25]. Although this is not a higher-order derivative 
problem it exhibits negative norm states. We show that using the reality conditions we can 
eliminate them. Furthermore, we evaluate the Green's function of the system using the Schwinger 
quantum action principle and the path integral method, finding that both of these procedures 
lead to the same result. In Section H] we combine the results of the previous sections in order 
to make a first attempt to extend the formalism to HOTD theories, successfully taking as an 
example the Pais-Uhlenbeck oscillator, which can be treated in a similar way as the complexified 
harmonic oscillator. In addition we perform the canonical analysis in an extended phase space 
in order to compare the cases of PU with equal and unequal frequencies. We observe that in the 
case of equal frequencies, the Lagrangian transformation induces a non canonical one which, at 
the quantum level, will lead to a non diagonalizable Hamiltonian. Finally, concluding remarks 
are presented in Section [5j Unless otherwise stated we use units in which h = 1. 



2 General considerations 

Let us consider a system whose Lagrangian contains up to the n-th order time derivative in 
a quadratic form. Written in normal coordinates this Lagrangian will have the form, 

L(x,i;,...,x(")) = ^ai(a;W)^ -y(x), i = l,...,n, (2.1) 

i 

with x^^^ = ^§ and the ctj G M being constants. We are assuming that under the transformation 
to normal coordinates we diagonalize the kinetic term and that the potential only depends on 
the coordinates. This Lagrangian yields an equation of motion of order (2n) 

Vf-i)'^— = 

k=0 
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Assuming x and hence L to be real leads to the well known problems of the higher-order 
derivative models. So, as a first step in our extension we consider the analytic continuation 
of the variable x to the complex plane. Next we assume that the Lagrangian (j2.ip plus a total 
time derivative will be real. This certainly does not alter the classical equations of motion, 
but modifies the Hamiltonian formulation by the addition of new terms in the momenta. Now, 
the essential points of our procedure are the following two i) we assume that there exists a 
nonlocal transformation from the Lagrangian (j2.ip plus a total time derivative to a second 
order real Lagrangian with n independent configuration variables i^*. Specifically we assume 
the complexified description can be related to an alternative one in terms of a real Lagrangian 

through 

L (x, X, . . . , X + I = . . . , C, r ) , (2.2) 

/ = /(x,...,xW), (2.3) 

r = r(x,...,x(")), i = l,...,n. (2.4) 

Here i = 1, . . . ,n are real coordinates and only their first order time derivatives enter 
in in a quadratic form. In general / will be complex, ii) The second point is that the 
transformation (j2.4p will impose reality conditions on the x^'^\ with /c = 0, . . . ,n. This implies, 
in particular, that terms in the potential of higher order than quadratic must be real functions 
of the variables, i.e., under the analytic continuation, we must consider only the real part of 
the interaction, for example powers like x^ should be replaced by (xx^)^. This is not necessary 
for the quadratic terms, since our transformation incorporates them properly. We must notice 
that strictly speaking, our system is different from the original one, since we have analytically 
continued the x variable. However, the reality conditions allow us to construct a properly well 
defined problem, with the same number of the original degrees of freedom together with the 
same classical equation for the original variable x. 



2.1 Canonical transformation and boundary conditions 

One of the interesting properties of our nonlocal transformation (12. 4p is that in the Hamiltonian- 
Ostrogradsky formalism [31], this transformation corresponds to a local canonical transforma- 
tion. To see this, we observe that in this formalism the corresponding phase space will be of 
dimension 2n, so the configuration space is extended io Qq = x^ Qi = x, . . . , Qn-i = x^^~^\ 
together with the momenta 

E (-^J . = 0,...,n-L (2.5) 

j=0 

Thus, the transformation (|2.4p only depends on the indicated phase space variables, 

r(x,...,x(")) =r(Qo,Qi,...,Qn-i,n„_i), (2.6) 

in which only the last momentum H„_i appears. We also have the corresponding transformation 
for the momenta, 

Pi = PiiQo, Qi, ■ ■ ■ , Qn-ijHo, . . . , n„_i), (2.7) 

given directly from the generating function (|2.3p . The transformations (|2.6p , (|2.7p will be canoni- 
cal in an extended phase space because the Lagrangians L and differ only by a total derivative 
that depends on the coordinates and momenta. We must remark that the transformations (|2.6p . 
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(|2.7p should be complex, because an analytic continuation of the original variables was required. 
However, these transformations define exactly the reality conditions, since we enforce that the 
2n variables will be real and, in consequence, also define the integration contour that 

must be used to quantize the system. 

To build the explicit form of the canonical transformation that we are using, we consider the 
equivalent of the expression (j2.2|) in Hamiltonian form. 



- n{Q, n) + ^^^^ = - H{t p). (2. 

From the comparison of independent terms we find 
df p df _ di^ 



These equations prescribe the 2n partial derivatives of the generating function f{Q,Il) and the 
Hamiltonians are equal due to the fact that the transformation is time independent. 

Now, to analyze the boundary conditions that are imposed by our complex canonical trans- 
formation, we make the variation of the left hand side of expression (j2.8p and we observe that 
the boundary condition is associated with the total derivative 

which dictates the appropriate combinations of the variations 6Q^ and Jllj to be fixed at the 
boundaries. Using the equations (j2.9|) we can show that the above expression is indeed ^ {PiS^^) , 
in terms of the variation of the So, the variables that we need to fix on the boundary are 
those functions of and Pi defined by the .^*'s. 



3 Reality conditions and the complex oscillator 

Let us start by reviewing the example of a modified harmonic oscillator subject to a complex 
canonical transformation [24:\ I25| . This has also been studied in a form that exploits the VT 
symmetry in the framework of the non-Hermitian Hamiltonian approach \23\ [321 [33l . The 
strategy is based on the observation that the complex Lagrangian 




where e is a real parameter, becomes the one corresponding to the harmonic oscillator after 
adding to it the total time derivative fc = In the Hamiltonian description the 

canonical momentum is 

dLc . . 
P = -jTr- = q- ieq, 
oq 

so p € C and should be quantized as a non-Hermitian operator. 

3.1 Hilbert space and reality conditions 

The quantum Hamiltonian becomes 



Hc = ^ + ^--/q' + ^{v,q}, 



(3.1) 
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where {p, q} = pq + qp so that a symmetric ordering is selected. Upon the canonical transfor- 



mation p = P — ieQ, q = Q, with P = -g-, the Hamiltonian ()3.ip becomes the standard one of 

the harmonic oscillator Huo ~ ^ ^ which is Hermitian whenever Q and P are. 
The coordinate representation 

qijj{q) = qipiq), p^{q) = -^^(9)1 

leads to a non-Hermitian form of the Hamiltonian (j3.ip with the usual scalar produc10. The 
corresponding Schrodinger equation becomes 

^"{q) - 2eq^'{q) - ((1 - e^)q^ + e - 2E) ^{q) = 0. 

Hence the eigenvalue problem for ()3.1|) can be related to that of the harmonic oscillator by using 

the wave functions of the latter = A'^e 2 Hn{q) and defining 

1 

il^niq) ■■= e ^ ^n, En = n+^, n = 0,1,2.... (3.2) 

A tedious but otherwise direct calculation shows that the eigenfunctions V'n do not have positive 
norm in the Hilbert space Tio = L^(M, dq). However, the classical role of fc as the generator of 
a canonical transformation motivates the introduction of a quantum transformation such that 
(see for example [Ml EB] for a general discussion) 

Hno = e-'"^ Hce'"^ . (3.3) 

2 

This non- unitary canonical transformation changes the measure dq to d/i = e~^^ dq. Now the 
Hilbert space TC = L^(M, d^) ensures (|3.2p have positive norm since 

(n|m)^ := j dqe'""^^ ip*'iprn = j dq^pl^Lpm =■ {n\m)g = d^n, 

where subscripts denote the appropriate measure for the correspondent Hilbert space. Moreover, 
the reality conditions 

q^ = q, p^ = p + 2ieq (3.4) 



are automatically implemented in 7i, where (13. ip is Hermitian 

It is illuminating to consider the Green's functions corresponding, respectively, to He 
and Huo- We recall they can be defined as 

Gh{q2,qi,E) := {q^\—l-^\q,) , (3.5) 
E — He 

GRo{q2,qi,E) := {q2\— ^ — \qi). (3.6) 

E — Huo 

Now on account of (j3.3p . we can rearrange (j3.5p to relate it with (|3.6p as follows 

Gh{q2,qi,E) = (g2i ^_^_l^^^^j gi) = (^^lor^^^^aki) 

= e+<1-'i^Guo{q2,qi,E), := , (3.7) 

where in the final line we use that q acts diagonally on the basis \q) and its dual. Thus the stan- 
dard harmonic oscillator Green's function, G}io{q2,qi, E), is obtained from that corresponding 

to the complex oscillator, Gq((72, qi,E), upon multiplying the former by e*^^ 2 2 ) , 



^Using p = — + ieg would yield the usual Hamiltonian for the harmonic oscillator which, however, requires 
a modified inner product to be Hermitian. 
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3.2 The Schwinger quantum action principle 

We can use Schwinger's quantum action principle to derive some dynamical properties of the 
system [37]. For instance, we can calculate the propagator {q2,t2\qi,ti) . We will work in the 
Heisenberg picture, where we adopt the notation 

{q2,t2\q{t2) = q2{q2,t2\, q{ti)\qi,ti) = qi\qi,ti) , (3.8) 
{pht2\pHt2) = P2{P2,i2\, p{ti)\pi,ti) = pi\pi,ti), (3.9) 

To achieve our purpose, we first calculate the total variation of {q2,t2\qi,ti) , given by 

Hq2,t2\qi,ti) = i{q2,t2\Gq^ - - H5T\qi,ti), 
where T = t2 — ti, and the generators Gq^ are given by 

Gq^ = p{t2)Sq2, Gq^ = p{ti)6qi. 

Starting from the Hamiltonian operator p.ip . we get the following Heisenberg operator equations 

q{t)=p{t)+ieq{t), 

p{t) =-{l-e^)q{t)-iep{t). 

Fixing the final and initial conditions at t2 and ti as q (t2) and q (ti), we obtain the solutions 

q{t) = A cos(t) + B sin(i), (3.10) 
p{t) = {B- ieA) cos(t) - {ieB + A) sin(t), (3.11) 

where the coefficients A and B are given by 

A = . ^ {q{ti) sin (i2) - 9(^2) sin {ti)) , 
sm (i j 

B = . {q{t2) cos (ti) - q{ti) cos {12)) . 
sm [T) 

In this way, we can express ^(^2) and p{ti) in terms of the initial and final conditions q{ti) 
and q{t2), getting 

p{h) = -^{-{iesm{T)+cos{T))q{ti)+qit2)), 
sm[I ) 

P{t2) = -4;-(-(iesin(T) - cos{T))q{t2) - q{h)). 
sm[i ) 

With this relations, one can also calculate the commutator 

[q{h),q{t2)]=ism{T). 
Using all of the above, we can write the Hamiltonian as 

H = ITT^p^, {Q\tl) + q\t2) -2cos{T)q{t2)q{h)) - 

2sm^(r) ^ 2 sm(T) 

Finally, from the equations (j3.8p . (|3.9p we get: 

S{q2,t2\qi,ti) = i{q2,t2\l | sin(r) - qi) - i£q2^Sq2 
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+ i sin(r) ~ ^2) + ieqi^ 5qi 

^ 'ql + ql - 2cos(r)g2gi) - \ST)\q,M 



= ^(92,t2|5|^-^^ {{q^ + qf) cos{T) - 2q2qi) 

Integrating we obtain 

{q2, h\q,,t,) = -^^J— exp {j-^^{i4 + 4) cos(r) - 2q,q,] 

xexp{^(gi-g?)}, (3.12) 



which is consistent with equation (j3.7p . Using the same procedure we now wish to obtain the 
quantity (^25^21^11^1)- Again, we begin by calculating its variation 

<J(P2i*2|pi,ii) = {p*2,t2\Gpi, -Gp^ -H6T\pi,ti), 

where now 

Gp* = -q{t2)Sp2, Gp^ = q{ti)5pi. 

We now need to fix the initial and final conditions at ti and ^2 as p{ti) and ^^(^2)- Thus, again 
we get (|3.1U|) and (|3.11|) as solutions to the Heisenberg equations, but now with coefficients A 
and B given by 

1 

A 



B 



((e2 + l)sin(T) - 2iecos(r)) 

X {(cos(t2) + iesni{t2))p{ti) - (cos(ti) - zesin(ti))p'''(t2)}, 
1 

" ((e2 + l)sin(r) - 2iecos(T)) 

x{ (sin(t2) - iecos{t2))p{ti) - (iecos(ti) + sin(ii))p'''(t2)}- 



In this way, we can express q{ti) and q{t2) in terms of the initial and final conditions p{ti] 
and ^^(^2)- The result is 

= ((.^ + l)sm(r) - 2..cos(T)) + ^^^^MTMH) - At.)], 

One can also calculate the commutator 

[p{ti),p\t2)\ = i ((e' + 1) sin(r) - 2iecos(r)) . 

Using all of the above, we find the following expression for H in terms of the operators p^(ti) 
and p(to) 

H ' 



2 ((e2 + 1) sin(r) - 2ie cos(r))' 
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X {(1 - e^){p\h)+p^^it2)) - 2((e2 + l)cos(T) + 2iesm{T))pHt2)p{ti)} 
i ((e2 + l) cos(T) + 2iesm(T)) 



2 ((e2 + 1) sin(T) - 2iecos(r)) " 
Since p^{t2) and p{ti) act according to equation (13. 9|) we get 

S{P2,t2\pi,tl) = i{p2,t2\' 



((e2 + l)sin(T) - 2iecos(r)) 
X |((cos(r) + iesin(T))p2 - Pi) + ((cos(T) + zesin(r))pi - p\) 5pi^ 

:{{l-e'){pl+pl')-2{{e^ + l)co^{T) 



1 



.2((e2 + l)sin(r) - 2iecos(r))^ 
^9- • * \ ' ((g' + l)cos(r) + 2iesin(r)) ^^^), 

= ^(P2,*2|5 2((e2 + l)sin(^)-2^.cos(^)) {(^"^^'^^ + ^^^^"^^^^ " '^^^4 



— ilnf — = I >|pi,ti). 

\ V(e2 + 1) sin(r) - 2ie cos(r) / J ' ' 

Integration produces 



(P2>*2|Pl,il 



V^(e2 + 1) sin(r) - 2ie cos(r) 

' i {icos{T)+iesm{T)){pf + pj) - 2p*^pi) 



X exp 



2 ((e2 + 1) sin(r) - 2iecos(r)) 
where C is a normalization constant. We note that 



lim {P2,t2\pi,ti) = CA/--^exp|--3-(pi - p*2f 

12^11 V 46 

Furthermore, in the hmit e — > 0, ah eigenvalues p cease to be complex, so in the double limit we 
must have 

lim(P2>*2 ti\pi,ti) = 5 {P2-Pi) ■ 

Fixing the constant C accordingly we obtain 

1 1 



{P2,t2\Pl,ti 



^ \/(e2 + l)sin(r) - 2iecos(r) 



[ i ((cos(T) + i6sin(T))(p^2 ^p2^ _ 2p*p^^, | 
^^^P|2 ((e2 + l)sin(r)-2iecos(r)) j' ^"^'"'^^ 

Next we compare the standard harmonic oscillator propagator with that obtained in equa- 
tion ()3.13p . The bracket corresponding to the change of basis is readily shown to be 
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which aUow us to relate the propagator expressed in terms of Hermitian variables, (P, Q), with 
that written in terms of non-Hermitian ones, {p, q), as follows 



We have explicitly verified the above expression. The completeness relation in the non-Hermitian 
description results 

Lp+i(p-p*)' 



d pKp^P*)\p){p*\ = 1> = 
3.3 Path integral approach 

Consider the action for the modified harmonic oscillator in terms of the real and imaginary 
parts of q and p, 

S = dt(^pq- (^^ + ^{1- e^)q^ + iep(^ ) = [ '^4(PR + {qR + Hi) 



- + ^(1 - e^) [qR + iqif + ie{pR + ipi){qR + ig/) j j • 

To compute the propagator using the path integral of this system we need to select an integration 
contour given by the reality conditions (|3.4p . that in terms of real and imaginary parts is 

qi = 0, pi = -eqR. 

These conditions are a pair of second class constraints in the framework of the Dirac's method 
of quantization [381 E^i so the path integral subjected to them is exactly the Senjanovic path 
integral [iQ]. In our case we obtain 



(g2, t2 151,*!)= j DpRDpiDqRDqi5 {qi)5 {pi + eqR) exp {iS) , 

since the determinant of the Poisson bracket of the constraints is one. Integrating over the 
imaginary parts using the delta functionals results in 

{q2,t2\qi,ti) = J DpRDqR(iy.Y>{i dt (^pR - ieqR) qR - + 

Here we see that the path integral is reduced to the usual path integral of the harmonic oscillator 
plus a term that contributes only to the classical action, so finally the amplitude is 

{q2M\qiM) = (^-^) ^ exp M + gi) COST - 2,1,^] + " , 

with T = t2 — ti. In this way we recover the result (I3.12P . 

Summarizing, we have provided a complex canonical transformation taking the complex har- 
monic oscillator (13. ip into the ordinary harmonic oscillator so that the non-Hermitian variables 
fulfill reality conditions in the Hilbert space 7i which are consistent with the canonical trans- 
formation. Let us observe that VT symmetry has not been invoked in this approach. The 
states (|3.2p possess positive definite norm in 7i and a unitary time evolution follows from the 
canonical transformation. 
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4 The modified PU oscillator as a second order theory 

Now we proceed with the HOTD model. Let us start with the PU oscihator Lagrangian 

Lpv = --x^ + ^ ^' x^ - -Y^x", (4.1) 
where we assume uoi > cj2- -^^pu is connected to 



1-2 Uj} r, 1-2 OJn o 

^5 = 7^6 -^ei + ^i2 -^ii (4.2) 



for real .^j, i = 1,2, by the following relations 

^pu + ^ = f = (^■^) 

S^i = i{ax + hx) , £,2 = cx + hx. (4.4) 

By choosing a, 6, c to be real we see that x is necessarily complex. To accomplish (14. 3p . 
the following values are obtained, for which the same sign should be used, 



1 



2 2 



V^'i - ^2 



(4.5) 



In other words, according to (|4.3|) . Lpu in (|4.1|) fails to be the real in (|4.2|) only by a total time 
derivative. Note that / will give rise to the corresponding canonical transformation once it is 
expressed in terms of phase space variables which we now derive. According to the Ostrogradsky 
method applied to (|4.ip we get 

Hx = [cvf + 0^2)^ + 'X , Z = X, liz = —X, 

SO that {x,!!^;} = 1 and {z,!!^} = 1. Clearly / = —zHz- The quantum Hamiltonian is 

^PU = -\fil - ^^^z' + + ^x\ (4.6) 

Using the afore mentioned transformation we have 

x = ihii + hi2, fix = iaPi + CP2, (4.7) 
z = ibPi + bP2, = icii + a^2, (4.8) 



where Pi, i = 1,2, are the canonical momenta conjugated to obtained from (|4.2|) . In terms of 
these Hermitian variables the Hamiltonian (|4.6p takes the form 

p2 , ,2 p2 , ,2 

H, = ^ + ^i! + f + f el. (4.9) 

So, starting from a Hamiltonian that is not bounded from below, extending the canonical vari- 
ables to the complex plane and adding a total time derivative a well defined Hamiltonian in 
terms of Pi and is obtained. Notice that in [16] it is acknowledged there exists a similarity 
transformation relating the original PU oscillator with a couple of independent harmonic oscil- 
lators similar to our case. Nevertheless a further transformation has to be supplied in [TU] to 
obtain the right frequencies for the oscillators. 
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4.1 The modified PU Lagrangian in the unequal frequency case 
in the extended phase space 

Instead of performing the canonical analysis starting from equation (j4.ip . we will consider (j4.3p 
which includes the additional time derivative. Namely we start from 

r ^ -2 , ^ f 2 , 2\ -2 ^ 2 2 2 I ■ (3) 

J^T = 2^ + 2 i^'i- ~^ ^ - 2^1^2X +xx^'. 

We first notice that Lx = -Z>t (x, i;, x, x^^^) , so it is higher order than Lpu- This also means 
that, according to Ostrogradsky's formalism, there will be a six-dimensional phase space with 
coordinates and momenta {pQ-,pi-,P2) given by 

= {ujI + ujI) ± + x^^\ (4.10) 
Pi = 0, (4.11) 
P2 = (4.12) 

according to the general definitions given in equation ()2.5p . Now we notice that (I4.1ip and (I4.12P 
are in fact primary constraints, which we write as 

<^i=m = 0, (4.13) 
(/>2=p2-i = 0. (4.14) 

The time evolution of the above constraints fixes the corresponding Lagrange multipliers making 
them second class constraints. In this way, the reduced phase space has dimension four. We will 
choose to express our Hamiltonian in terms of (x, x,po5P2)- The resulting canonical Hamiltonian 
is given by 

^^c = Pi)X + pix + p2X^^^ - Lt = PqP2 - i^x^ - ^ (^^1 + ^DpI + ^^i^W- 

Next we construct the canonical structure in the Dirac formalism, accounting for the second 
class constraints. To do so, we use the Dirac brackets 

{A, B}* = {A, B} - {A, C"^ {0^, B}, a, /? = 1, 2, 
where is a square antisymmetric matrix such that 

Upon using equations (j4.13p and (|4.14p we obtain 

The resulting Dirac brackets are 

{x,po}* = l, {x,p2}* = l, (4.15) 
and zero for any other combination. Next we rewrite the Lagrangian in the canonical formalism 

Lt = xpo + xpi + x^^^p2 - He 
and make the constraints (j4.13p . (j4.14p strong obtaining 

Lt = xpo + x(^)p2 - Hq. (4.16) 
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In order to make contact with the classical version of the Hamiltonian (|4.6p we introduce the 
following canonically related new variables, leaving x unchanged, 

^x=Po, z=p2, Uz = -x. 
Using equation (j4.15p . we obtain as the only non-zero Dirac brackets 

{x,n,}* = i, {z,u,}* = i. 

In this way, we can consider the variables {x, z,Ilx,Ilz} as a set of canonical variables in the 
Dirac formalism. Introducing the new variables in equation ()4.16p we obtain 

Lt = xU^ - ti-^z - He, 

where Hq corresponds to the classical counterpart of i/pu given in equation ()4.6p . The above 
equation guarantees that -ffpu is the correct Hamiltonian in the reduced physical space and 
imply that the correct boundary conditions at the end points are given by fixing x and II^ in 
regard to the variational problem. 



4.2 The Pais— Uhlenbeck transformation in matrix form 

The complex canonical transformation between the real variables {^1,^2, Pi, P2} and the complex 
ones {x, z, n^:, n^}, is given in matrix form by 

X = M^, C = M~^X, 

where X, ^ and M stand for 



X 



( X \ 

z 



( ?! \ 

6 
V ) 



M 



( ih b \ 

i6 6 

c 

\ic a / 



Now, the matrix M satisfies 

M^nM = n, 

with 



/ 1 \ 
001 
-1000 
V -100/ 



(4.17) 



where only the relation 6(c — a) = 1 was needed to prove equation (j4.17p . This means that the 
matrix M belongs to the complex symplectic group Sp(2n,C), with n = 2. Since Sp(2?i,C) is 
a 2n(2n-|- 1) parameter group, in our case M should have 20 parameters in general. One notices 
that 



det(M) = (6(c- a))^ 



1, 



which verifies the unimodular character of the symplectic matrix. 
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4.3 Hilbert space 

We have related the Hamiltonians (j4.6p and (j4.9p by means of the complex canonical transfor- 
mations (j4.7p . (j4.8p and their inverse 

^1 = iax — ihliz, S,2 = cx — bUz, 

P^-^ = —icz + ibilx, -P^2 — + bllx, 

where the coefficients a, b and c are given in ()4.5p . We notice that in this case the following 
relations are satisfied: 

ac = cvfoj'^b'^, b {c — a) = 1. 

At the quantum mechanical level, we fix the reality conditions by requiring that the operators 
{Cii 5 -Pgi) -^^2} Hermitian. This in turn implies that: 

x'^ = b{a + c)x -2b'^Uz, = -b{a + c)z + 2b^Ux, (4.18) 

Ul = 2acx -b{a + c)Uz, Ul = -2acz + b{a + c)Ux. (4.19) 

We will quantize the system in the basis \x, Hz), so we propose the following operator realizations 

X = X, n, = n,, 



and 



llx = —ll^, Z = l- 



dx dliz 

Now, the inner product written in the basis Ix,!!^) can be expressed in general as 



{(k, ij) = j (fxcPiiz fi (x, n,) r (x, n,) ^ (x, n,) , 



where 

(fx = dxfidxj, cfHz = dllzRdllzi- 
Because of conditions ()4.18p and (j4.19p . we should have in particular: 

(x(/-, V') = (</>, (6 (a + c) X - 2b'^Uz) ip), (4.20) 

and 

[Hzct), ij) = ((/., (2acx -b{a + c)Uz)il^). (4.21) 



or 



J (fx(fUzHX*(p*i; = j (fxcPYlziJ.(i)*{b{a + c)x-2b^Ilz}i}, 
J (fx(fUzfin*z(j)*ij = j (fx(fUzH(j)*{2acx-b{a + c) H^} V'- 



Separating x and Hz in their real and imaginary parts inside the integral, and matching both 
equations, we conclude from (I4.20p and ()4.2ip that: 

axR = bUzR, (4.22) 
cxj = bUzi, (4.23) 

which can be always accomplished if 

fi (x, n^) = C6 {axR - blizR) 5 {cxi - blizi) , 

where C is a constant to be determined by some normalization condition. 
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4.4 The Schrodinger equation 

Let us consider the stationary Schodinger equation 

HMij{x,n,) = E^{x,u,), 

or more explicitly, using the operator realizations given above 



1 o , 1/9 ON ^'^^h 1,9,, ^'^^h 



Using the relation between {^1,^2} and {x,!!^}, applying the chain rule 

d sr-^ dS,i d . d d 
dx ^ dx d£i dfi 8^2 ' 

d _ sr^ _ _ _^ _ , d 

and substituting we arrive at 

with regular solutions: 

V'm,n (6,6) = exp (wi^i + ^2^1) I -f^m -f^n • 

Let us now consider the ground state, m = ?i = 

V'o.o (6,6) = (^iCi +^2^2) 

or written in terms of x, 11^ 

V'o.o (a;, n^) = exp I ^ ( (cjia^ - uJ2<?^ + 2b {uj2C - ujia) II^x + {oji - LV2) H^) | . 
Separating x and II^ in their real and imaginary parts we get 

V'co (2;,n^) = exp|i(^(wia2 -W2C^) (xR-xj) + 2b {uJ2C - uJia) {U^rxr - U^ixi) 

+ b^ (oJi - UJ2) {uIr - nlj)) I exp {i<P (x, n,)} , (4.24) 

where (j)(x,Ilz) is a real function of x/j, xj, HzR and H^/ and therefore is the imaginary part 
of ipQfi written as a phase factor. It is not easy to see from (14.24p if the ground state function 
vanishes at |x| — > ±00, III^I — > ±00. However, substituting the reality conditions (|4.22p . (j4.23p 
we have 

I r -n \ I ^2 {ujI - ^i) 2 1 / ^1 (^1 - ^2) 2I 
^0,0 [z, n^) = exp <^ '-xji ^ exp <^ Xj '> , 

where we verify that the ground state function indeed vanishes at infinity. 

Finally, in the Appendix we compute the full propagator of the PU oscillator. 
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4.5 The modified PU Lagrangian in the equal frequency case 
as a third order system 

We now wish to follow the same procedure as in Subsection 14.11 for the restriction of the PU 
Lagrangian to the equal frequency case. Recapitulating, we begin from the corresponding equal 
frequency PU Lagrangian Lpu plus a total time derivative 

i, = ir„ + | = -l.^ + .^*^-i.V + f (4.25) 

The equations of motion arising from Lpu are 

= x(^) +2a;2x(2) +iv'^x, 

which are not affected by the total time derivative. 
At the Lagrangian level the choice 



f = - {x-u x)x 
allows to transform (j4.25p into the two oscillators system 

+ ft=k = \il + \il - ifi + ei) - 

by means of the nonlocal transformation 

£,1 = i {ax + bx) , £,2 = cx + dx. 
together with the conditions 

^^2 = 5^2^ 1)2^^2^ cd-ab= 1/2, (?-c? = ij?. 

For simplicity we can take d = b, which gives us 

a = a;^6--^-, c = u?b+-^. 

4b 46 

Then, in this case we have a free parameter in the transformation. Now, the total Lagrangian 
turns out to be 

f -^2-2 142,-^ fS"! • ^ 2 ■■ 
Lt = —UJ X LO X H X^ 'x to XX, 

2 2 2 2 

where x*-^-* = 'x. According to Ostrogradsky, the canonical momenta conjugated to the coordi- 
nates (x, X, x) are 

Po = -Lo'^x + x'-^\ (4.26) 

pi = -^{uj^x + x) , (4.27) 

P2 = h, (4.28) 
respectively. Again, (j4.27p and (|4.28|) are primary constraints, which can be put as 

(Ai =Pi + ^(w^2; + x) = 0, ,/,2 = P2 - = 0. (4.29) 
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As in the previous case, the time evolution of them fixes the corresponding Lagrange multiphers, 
thus turning (|4.29|) into a set of second class constraints which reduce the dimension of the phase 
space from six to four. We choose to express our Hamiltonian in terms of variables {x, x,po,p2}, 
eliminating x, pi from the constraints. The result is 



^^C = POX + PlX + P2X 



(3) 



Lt = 2p2{po + UJ2P2) - 4w^p| + ^oj^x"^ - ]^x^- 



A similar procedure as the one described in Subsection 14. II yields the non-zero Dirac brackets 



{x,po}* = 1, 



1 



-UJ 



{x,P2y 



1 

2' 



in the phase space {x, x^pQ,p2} ■ In order to make contact with the original equal frequency PU 
problem, let us introduce the new variables x, z, lix, 



X = X, 



2P2, 



tlx=PQ + U?P2, 



and proceed to calculate their Dirac brackets. The only non-zero results in the phase space 
{x^z^Ilx^^z] are 

{x,n4* = i, {z,n4* = i, 

The above indicates that we can consider the set of variables {x, 5, ITa;, IT^} as canonical variables 
in the Dirac formalism. Substituting in Hq we get: 



He 



■ uP'z^ -|- Haji H — lo'^x? 



(4.30) 



which corresponds to the original Pais-Uhlenbeck Hamiltonian (j4.6p in the equal-frequency case. 
As pointed out in references [El [21], (I4.30p is not diagonalizable and hence cannot be related 
by a similarity transformation to the two oscillator system. 

However, the Pais-Uhlenbeck Hamiltonian in the equal frequencies case can be related to the 
Hamiltonian of two uncoupled oscillators 



1 



1 



-Ijj 



through a transformation which is neither canonical nor of the similarity type, given by 



/ X 






z 














J 





ib 










id 


d 





i{dw^-h) 


{diu^ + 




bu^ 








/ 6 \ 








Pi 




\ P2 J 



or, conversely 



/ ^1 












Pi 






V P2 


) 





ibu 



ijj 



\ 



■4 + 62 











+ 



+ id 








62 \ 



/ X \ 

z 

V n. J 



Even though we have not explored the quantization of the system in this case, we expect that 
the non existence of a similarity transformation between the PU oscillator and the two-oscillator 
system would be a signal of the non-diagonalizable property of the Hamiltonian related to the 
appearance of Jordan blocks, as discussed in [HI [TTl [21]. 
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5 Discussion 

In this work we have proposed that a complex canonical transformation applied to a class of 
HOTD models with Lagrangians modified by a total time derivative can solve the drawbacks 
of such theories, namely negative norm states or ghosts, unbounded below Hamiltonian and 
non-unitary time evolution, keeping unmodified the original classical equations of motion. Such 
a transformation requires the original canonical variables to be analytically continued and fur- 
ther subjected to appropriate reality conditions so that the final canonical set matches the 
original number of real degrees of freedom. In particular while the final canonical variables ^* 
are Hermitian the original ones x^^^ fulfill some reality conditions dictated by the canonical 
transformation among such variables. The inner product in the original Hilbert space is defined 
in such way to reproduce the reality conditions. This idea has been applied previously to non- 
perturbative canonical quantum general relativity [24^ [25] and the inclusion of reality conditions 
was mandatory to recover real gravity. Thus we have successfully provided a novel extension 
of such idea to the case of the PU oscillator as an example of HOTD models. An important 
point that we must clarify is that the analytic extension of the variables allow us to have a well 
defined transformation between the HOTD model and a second order theory. The interesting 
result is that both theories have the same degrees of freedom and are related by a complex 
transformation. So, our procedure to extend the original problem to the complex plane gives us 
the appropriate frame to find the mapping between both theories. The success of the proposed 
approach in the PU case motivates the study of its extension to other HOTD models, but this 
is beyond the scope of the present work. 

To illustrate the procedure we have first reviewed the case of a complexified model that 
reduces to an harmonic oscillator. Although not of the HOTD type, this example exhibits 
negative norm states which are neatly eliminated using a complex canonical transformation. In 
this case we determined also the corresponding propagator and contrasted it with the usual one 
of the harmonic oscillator. To do so we adopted two different techniques: Schwinger's quantum 
action principle and the path integral approach, both of which led to the same result. Let us 
notice also that the path integral version of complex canonical variables for gravity has been 
studied [i2] . 

Next we studied the PU oscillator ()4.ip which is properly a HOTD model. After complexi- 
fying the original variables we constructed the complex canonical transformation (j4.7p . (j4.8p 
which leads to a system of two decoupled harmonic oscillators ()4.9p with just the frequencies u>i 
and L02 appearing in (j4.ip . The propagator of the model containing non-Hermitian variables was 
determined by adopting again Schwinger's quantum action principle. The equal frequency PU 
oscillator has been also described in the classical case. We show that there are some differences 
with the unequal frequency situation. First, we have here a free parameter in the transformation, 
perhaps a resemblance of the Barbero-Immirzi parameter in loop quantum gravity |29l I30j . 
Furthermore, in the reduced phase space there does not exist a similarity transformation between 
the PU oscillator with equal frequency and the second order oscillators. However, we found that 
this system is also mapped to a couple of harmonic oscillators, using a non-similarity type 
transformation. 

Our approach to the PU oscillator is an alternative to that of [16] based on non-Hermitian but 
VT symmetric Hamiltonians. We find ours has the following properties: (1) It is based on non- 
Hermitian variables subjected to specific and well defined reality conditions arising from (j4.7p . 
(14. Sp that make the higher-order Hamiltonian Hermitian in the appropriate Hilbert space with 
measure ()A.3p . (2) No symmetry is required. (3) It works for at least the PU oscillator. (4) If 
we want to consider additional interaction terms as an anharmonic contribution, our proposal 
still works provided that we add the prescription {xx^ 

)2 = ^4(^2 _^ ^2)2^ rpj^-g interaction 
term has the correct signs and the Hamiltonian is still bounded from below. (5) We must 
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emphasize that the complexification of the original problem together with the reality conditions 
are a key ingredient in our approach which constitute the prescription that define our problem. 

It is remarkable that this approach can be also applied to the case of a quadratic higher-order 
derivative free scalar field model. To see this notice that the following Lagrangian density 

is related to 

= -iviD^i - ^V'? - ^V'2n^2 - ^Vi, (5.1) 

by 

ipi = i{a(j) + fen</'), il)2 = ccf) + ftD'/', 

with o, 6, c having the same form as in the PU model ()4.5p except for the replacing uji rrii, 
i = 1,2. Clearly (jS.ip is ghost-free. Our approach seems promising even in the interacting case 
including a term of the type {(t>4>*)'^ = — , ^ (V^? + V'i)^ ElE], however all these generalizations 
deserve further study. 

In closing we would like to mention some directions for further study along the lines proposed 
in the present work. Other physical mechanical models including at least third order derivatives 
might help to show the strength of the canonical approach here pursued. Also a crucial problem 
to come to terms with is field theory. The success in the simplest non interacting scalar field 
briefly described above requires full implementation at the quantum level. Finally, as it was 
mentioned in the introduction, higher order gravity has been one of the motivations to the 
present approach, which of course is a paramount challenge. Also it would be rather interesting 
to look for HOTD gauge field theories to test some aspects of the gravitational case. 



A Schwinger's quantum action principle and the PU oscillator 

In this appendix we compute the propagator of the PU oscillator, a result already obtained 
in [43]. To do this, we apply Schwinger's quantum action principle to calculate (rrg, 11*2, t2|3;i, 
n^i,ii). However, in this case we need to be more careful since our variables are complex. We 
will work in the Heisenberg picture, where the involved operators act as: 

x{ti)\xi,ii;,i,ti) = xi|xi,n2i,ti), i\z{ti)\xi,iizi,ti) = n^i|2;i,n^i,ti), 
(4,n:2,t2kt(t2) = {xi^ni^M^l, {xl,K2M'^\{i2) = {xi,i\':2MK2- 

As usual, we first calculate the total variation of the quantity (x2, 11*2, t2|3;i) n^i, ti), given by: 
5{xl,K2Axi,KiM) (A.l) 

= i(x^,n';2,t2|{4(t2)fe^ - z\t2)5K2 - n,(ti)fei + - ^5r}|xi,n,i,ii), 

where T = t2 — ti. Now, we must express {z, Ila;, z\t),Yi\{t)} in terms of {x, Yiz:x\t),Ii\{t)} . 
We solve the Heisenberg equations for the operators x, z, H^; and H^. This equations are: 

±{t) = z{t), 

^xit) = -UjluJ2x{t), 

z{t) = -n,(t). 
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U,{t) = {iof + u;l)z{t)-U^{t). 
Uncoupling, we get the following equation for x{t): 
+ {loI + ujj) x{t) + ujlujlx{t) = 0, 
with general solution: 

x{t) = C7ie^^i* + Cse-^"^!* + Cge^'^^i + c^e"''^'^ 
while for the other operators we get: 

z{t) = iuJi (Cie^"^!* - C2e"*^i*) + W2 (Cge^'^^t _ de-'"^^*) , 

We fix the constants Cj's in terms of initial conditions at time ti: x{ti), 112(^1) and final 
conditions at time t2- x^{t2) and Ill{t2)- To obtain such conditions, we use the reality properties 
(|4.18p . (j4.19p . valid at all times. In this way, we arrive to the following system of linear equations: 

x{ti) = C7ie*'^i*i + Cse"*'^!*^ + Cae'"^^^' + de-''^^^^ 

x^t2) = - (Cie^'^i*2 ^ C2e-''^'^^) + (C73e*^2t2 ^ c^e'"'^^^^) , 

Solving for the C's and substituting, we get the following expressions for z{ti), Ilx{ti), z^{t2) 
and ni(t2): 

z{ti) 



{loi — (JJ2) sin(a;iT) sin(u;2T) 

X {{—ujfuJ2 sin(a;iT) cos(a;27') + "^i"^! cos(lo'iT) sin(a;2T)}x(ti) 
+ {lo2 sin(wir) cos(t(J2T) — wi cos(u;iT) sm{uj2T)}Ilz{ti) 

+ {ujIuj2 sin(6JiT) + ujiujI sin(a;2T)}x^(t2) + {-UJ2 sin(6JiT) - uJi sin(i^2r)}nj(i2)}, 
1 

{ujf — sin(a;iT) sm{i02T) 

X {{—ujfuj2 sin(ci;iT) cos(a;2^) + ^1^2 cos(uJiT) sin(u;2r)}2;(ti) 

+ {ujIuJ2 sin(a;iT) cos(a;2T) — ti^iwl cos(a;iT) sin{uJ2T)}Ilz{ti) 

+ {a;fa;2 sin(cJiT) + UJ1UJ2 sin(a;2T)}x^(t2) 

+ {-ujIuJ2 sin(a;it2) - ^i^l sin(a;2T)}n|(t2)}, 
1 

{ujI — 0J2) sin(a;iT) sin(u;2T') 
+ {ijJ2 sin(6<JiT) + uji sin(a;2r)}n2(ti) 

+ {u>fu>2 sin(LiJiT) cos(a;2r) — u>iu>2 cos{loiT) sin(a;27')}a:^^ (^2) 
+ {-UJ2 sin(i^ir) cos(w2T) + uJi cos(wiT) sm{uj2T)}Ul{t2)} , 

nU*2) = 1-^ 2^ ■ / w rT.A {-^i^2sm{uJiT) - uJiUJ^sm{uJ2T)}x{ti] 

[ujf — LO2) sm(u;iJ j sm(w2-' j 

+ {uftjj2 sin(u;iT) + witt^l s'^'^{^2T)}Ilz{ti) 



-Z^(*2) = -;—2 2\ ■ I rr\ ;7;v{{~'^l'^2 sin(wiT) — UJ\Lo\'SVCL{uJ2Ty\x(t\ 
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+ {ujfi02 sin(a;iT) cos{uJ2T) — Wiwl cos(a;ir) sm{uj2T)}x\t2) 
+ {—uj\oj2 sm(wiT) cos{uJ2T) + iOiio\ cos(t(JiT) sm.{L02T)}Ti\{t2)} ■ 

Using these results in the Hamiltonian (j4.6p and in the total variation (jA.ip . we get finally 



(x2,n*2,t2|2;i,n^i,ti) = — . . , — ^exp 



sin(a;iT) sin(a;2T) 1 [lo\ — Wg) sin(u;ir) sin(a;2r) 
X < {u)\uj2 sin(Lo'iT) cos{uJ2T) — ujiuj2 cos(i^iT) sin(Lo'2r)} "j ^ "I — 1~ r 

f n* 2 

+ {u!2 sin(Lo'iT) 005(1^2^) — ^1 cos(a;iT) sin(a;2r)} < -\ ^ 

+ {u;^a;2 sin(a;iT) + wia;| s\n{uJ2T)}{xlUzi +11* 2X1} 

+ {—uj'fuJ2 sin(u;iT) cos(u;2T') + cos(wiT) sin(LL'2T)}{x2n*2 + xill^ 1} 

+ {—ujfuJ2 sin(u;iT) — wiwl sin(a;2T)}{x2Xi} 

+ {-L02 sin(ix;ir) - uji sin(u;2r)}{n* 1 1 ' 

or, for short: 



(x^ , n';2 , t2 1 XI , n, 1 , ti ) = g exp I - 1 F(x*n, 1 + xin: 2 ) + en, in: 2 

+ j(x^n:2 + xin,i) + K(^Si±^) +Mx,x*2 + n(^^^±^^^Y (a.2) 

with D, F, G, J, K, M, N , Q, being the following functions of T: 

D = {uj\ — UJ2) sin(u;ir) sin(a;2T), 
F = ujIuj2 sin(ix'ir) + a;ia;| sin(a;2T), 
G = —UJ2 sin(a;iT) — oji sin(a;2r), 

J = —Lo\uj2 sin(Lo'iT) cos{uJ2T) + uJiUj\ sin(tJ2^) cos(a;iT), 
K = UJ2 sin(LiJiT) 005(1^27") —oJi sin(tiJ2T) cos(u;ir), 
M = —ijj\oJ2sm{uJiT) — wi^l sin(t(J2T), 

= uj\uj2 sin(u;iT) cos(LiJ2r) — 1^11^2 sin(ti;2r) cos(a;ir), 



Q 



1 



sin(t(JiT) sni{uj2T) 



Just as in our first example the PU propagator ()A.2p can be related to the one corresponding 
to the Hamiltonian (14. 9p by using the change of basis 



(Pi, P2|x,n^) = exp [(ax - 6n^)Pi + {-icx + ibIl^)P2] ■ 
The basis [xjll^) is complete with the measure 
dxRdxidU^RdU^i 

dfip^ = -— ^ o[bnzR - axR)d[bUzi - cxi), (A. 3) 

where xr, xj, II^r, HzI are real. 
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